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On July 4, 1997, the Mars Pathfinder (MPF) mission successfully landed on Mars. 
The entry, descent, and landing (EDL) scenario employed the use of a Disk-Gap-Band 
parachute design to decelerate the Lander. Flight reconstruction of the entry using MPF 
flight accelerometer data revealed that the MPF parachute decelerated faster than 
predicted. In the summer of 2003, the Mars Exploration Rover (MER) mission will send 
two Landers to the surface of Mars arriving in January 2004. The MER mission utilizes a 
similar EDL scenario and parachute design as that employed by MPF. As a result, 
characterizing the degree of underperformance of the MPF parachute system is critical 
for the MER EDL trajectory design. This paper provides an overview of the methodology 
utilized to estimate the MPF parachute drag coefficient as experienced on Mars. 
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Background 

On July 4, 1997, the Mars Pathfinder (MPF) 
mission successfully landed on Mars. The entry, 
descent, and landing (EDL) scenario employed by the 
mission is shown in Fig. 1. Mars Pathfinder utilized 
numerous methods for decelerating the Lander during 
the entry. Early in the descent, the Lander was 
decelerated solely by the use of aerodynamic drag 
generated by the capsule forebody. At approximately 
8.6 km altitude, a parachute was deployed to further 
decelerate the Lander during the descent. 
Approximately at 15m above ground level (AGL), retro- 
rockets were fired to remove any remaining velocity at 
which point the Lander was released. Prior to Lander 
release, airbags were inflated around the Lander to 
cushion it during ground impact. References 1 gives an 



overview of the Mars Pathfinder mission, as well as, a 
description of the EDL events. 


1) Direct Entry from Hyperbolic Approach 

2) Cruise Stage Separation: E- 30 minutes 

3) Atmospheric Entry: -125 km altitude 

4) Parachute Deploy: -8.6 km AGL, -E+ 173 s 

5) Heatshield Jettison: 20 s after chute deploy 

6) Bridle Descent: 1 0 s after heatshield jettison 

7) Radar Acquisition of Ground: -1 .5 km AGL 

8) Airbag Inflate: -4 s prior to retrorocket ignition 

9) Rocket Ignition: -90 m AGL 

10) Bridle Cut: -15 m AGL, m/s vertical velocity 

11) First Contact w/Ground: -E+ 305 s 




Figure 1. MPF Entry, Descent, and Fanding Sequence 
of Events. 


reexamined. In addition, an overview of the 
methodology utilized in this reexamination of the MPF 
parachute Cd estimate is provided. 

Approach 

At the terminal velocity condition, the sum of the 
vehicle drag force (F D ) and the buoyancy force (F B ) 
equal the weight of the vehicle (W) 



The Pathfinder mission utilized a Disk-Gap-Band 
(DGB) parachute design to provide deceleration during 
the terminal phase of the entry. The DGB parachute 
used by MPF was a derivative of that flown on the 
Viking missions. 2 Reference 3 provides a description of 
the design and development of the MPF parachute 
system. 

Flight reconstruction of the entry using MPF flight 
accelerometer data revealed that Pathfinder decelerated 
faster than predicted based on the estimated value of the 
MPF parachute C D of 0.50; a value which was 
determined from low altitude Earth flight tests and wind 
tunnel data during the development of the MPF 
parachute (see Ref. 3). An explanation of this 
underperformance of the MPF parachute system from 
that which was predicted is still not known. 


W = mg (1) 

F B = pg Vol (2) 

F d = P y2 (CD Par Ap ar + C Db/s A b/s + C DLan A Lan )/2 (3) 

where 

m = suspended mass, 
g = gravitational acceleration 
p = atmospheric density, 

Vol = total vehicle volume 
v = terminal velocity, 

C Db/s = Backshell drag coefficient, 

A b/s = Backshell reference area 
C D = Lander drag coefficient, 

Lan 

A Lan = Lander Reference area 
Co p = Parachute drag coefficient, 

Apa, = Parachute reference area 


In the summer of 2003, the Mars Exploration Rover 
(MER) mission will launch two Landers that arrive at 
Mars in January 2004. The MER mission utilizes a 
similar EDL scenario and parachute design as that 
employed by MPF. Therefore, characterizing the degree 
of underperformance of the MPF parachute system is 
critical for the MER mission, since the parachute drag 
coefficient (Cd) is a driver for the MER entry, descent, 
and landing trajectory design. As a result, a good 
estimate of the MPF parachute C D as experienced at 
Mars is needed, along with a measure of its uncertainty 
in order to robustly design the MER EDL phase of the 
mission. This investigation was initiated to determine 
the best possible flight reconstruction estimate for the 
MPF parachute Cd as exhibited on Mars. All data 
sources required for calculation of the parachute Cd are 


Equating the forces and substituting for the 
parameter yields the equation for predicting the 
parachute C D 

CD Par = 2(mg - pgVol)/( pv~A Par ) - C Db/s A B / S /A P ar 

- CD La] A Lan /. Apar (4) 

Note, the velocity and density are the largest 
contributors for determining Cn p , since the velocity 
term appears to the second power in the equation, and 
any uncertainty in density will directly relate to an 
uncertainty in C Dp because density and Cd are coupled 
together for determining drag. Hence, an accurate 
estimate for both parameters is necessary. 
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All parameters required for calculation of the 
parachute C D are reexamined: MPF altimeter flight data, 
atmospheric properties, Backshell and Lander Cd 
values, parachute area, and total vehicle volume. Values 
for the parameters are updated where appropriate. 
Specifically, outlined is the rationale for the 
improvement in the estimates for the value of the 
terminal velocity and atmospheric properties, which are 
the largest contributors to the parachute C D estimate. 
Note, conservative assumptions are utilized where 
appropriate. 

To estimate the MPF parachute Cd value, this 
investigation employs the use of a Monte Carlo 
technique. In so doing, uncertainties in all the 
parameters can be included in the analysis. As a result, a 
statistical range on the MPF parachute C D can be 
defined. With this approach, a more rigorous 
methodology is employed which provides for the best 
possible flight reconstruction estimate for the MPF 
parachute Cd- 

Atmospheric Density Determination 

Accurate density estimation is necessary for 
determining the parachute C D , since the two parameters 
are coupled. However, neither density, pressure, nor 
temperature were measured during the MPF parachute 
descent. 

Density can be derived through temperature and 
pressure measurements. Accurate MPF temperature and 
pressure measurements are available at the surface 
(altitude of 1 m). However, measurements are required 
above the surface near terminal velocity conditions. 
Therefore, surface pressure and temperature must be 
extrapolated to the altitude where C D is derived. 
Pressure can be extrapolated accurately, however, 
temperature must be derived from models or other 
measurements. The corresponding density can then be 
calculated by extrapolating the pressure and 
temperature, which were measured at 1 m. The 
temperature increases rapidly with altitude from the 
surface over the lower few hundred meters. However, 
considerable uncertainty exists for extrapolating the 
temperature from the surface to higher altitudes due to 
the uncertainty in the width of the temperature inversion 
layer. 

Therefore, current atmospheric property estimates 
are supplemented by Mars Global Surveyor (MGS) 
radio-occultation data. 4 See Appendix for a detailed 
description on the rationale. Current best estimate of 
MPF parachute Cd is performed at an altitude of 1 km. 
The temperature at 1 km is taken from profiles derived 
from MGS radio-occultation measurements. High 


vertical resolution profiles from the surface to 30 km are 
obtained near the Pathfinder season, local time, and 
location, but for the following Mars year. The profiles 
show a maximum temperature of 220 K at 1 km with 
little day-to-day variation. This estimate is based on the 
assumption that year-to-year temperature variability at 1 
km is small. 

The year-to-year temperature variability at 1 km is 
well within +9 K (3-0) error assumed. Multi-year radio- 
occultation measurements are not available for the MPF 
conditions. However, MGS Thermal Emission 
Spectrometer (TES) data show little year-year 
variability in mean 0-10 km temperature at other 
seasons. Multi-year TES coverage of the MPF season 
are not yet available, although. Viking and MPF 
Landers demonstrate that surface pressure varies little 
from year-to-year. However, due to the lack of accurate 
temperature data, a uniform distribution for the 
temperature variation is utilized in the analysis for 
conservatism, rather than a normal distribution which is 
more consistent with the observed temperature data. 

The assumptions outlined in this investigation 
regarding the atmospheric properties were 
independently reviewed. The approach taken is 
consistent with known modeling constraints and 
uncertainties. Overall, the approach developed is the 
best that can be constructed with the available data to 
define the atmospheric properties during the MPF 
parachute descent. 

Terminal Velocity Determination 

The MPF altimeter flight data was reexamined 
using higher order filtering techniques beyond that 
performed in the early reconstruction efforts 5 . The 
current investigation performs two independent filtering 
techniques in an effort to smooth the MPF altimeter 
flight data to obtain a more accurate determination of 
the terminal velocity. A Butterworth filter and a more 
sophisticated Kalman filter were utilized on the flight 
data. 

The Kalman filter (actually a backward-in-time 
filter plus a forward-in-time smoother) was designed to 
process data simultaneously from the MPF 3-axis 
accelerometer and the radar altimeter data. Nine 
parameters were estimated in the filter. These were: 
height above the ground, horizontal distance from the 
touchdown point, vertical and horizontal speed relative 
to the ground, vertical and horizontal non-conservative 
(i.e. drag) acceleration components and the time 
derivatives of these accelerations, and the angle between 
the MPF Lander base petal and the vertical direction. 
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The filter was first run backwards in time over the 
accelerometer and radar data set, from a point in time 
just prior to firing of the MPF retrorockets, to the time 
of the first radar measurement (a time span of 
approximately 20 seconds). The filter state vector was 
initialized by integrating the accelerometer data 
backwards in time, assuming a constant gravitational 
acceleration. The backward-in-time sequence of events 
included the initial impact with the ground, the free-fall 
of the Lander, the severing of the MPF bridle 
connecting the Lander with Backshell and parachute, 
and the firing of the retrorockets. A large covariance 
was used for the backward-in-time filter pass. 

Subsequently, the forward-in-time filter was 
initialized with the state vector of the backward-in-time 
filter, taken at the first radar measurement time. The 
state vector histories verses time of the forward and 
backward filters were averaged to produce the smoothed 
estimated state. The measurement residuals computed 
with the smoothed state had essentially zero bias, and 
RMS values of 1.47 meters for the altimeter 
measurements, and approximately 0.05 Earth g’s per 
axis for the accelerometer measurements. One 
component of the smoothed estimated state, the vertical 
velocity, is shown in Figs. 2 and 3. 

The two different filtering techniques show a good 
agreement in the mean value for the terminal velocity as 
seen in Figs. 2 and 3. In addition, as a by-product of the 
Kalman filtering technique, a statistical estimate on the 
uncertainty in the terminal velocity is obtained of +1.8 
m/s (3-0). This uncertainty estimate can be directly used 
in the Monte Carlo approach for specifying a dispersion 
on the terminal velocity. 



Figure 2. MPF Terminal Velocity Estimate using 
Filtering Techniques 



Figure 3. MPF Terminal Velocity Estimate using 
Filtering Techniques (zoomed in). 


Other Parameters 

The MPF parachute area (of the flight unit) was a 
quantity that was not explicitly measured. 
Unfortunately, there was no requirement to control the 
actual MPF canopy area. As a result, its exact value is 
unknown. Since the parachute area is a quantity that 
directly affects the value of the parachute C D , an 
accurate knowledge is required to calculate the value of 
the MPF parachute C D . 

To improve the knowledge of the MPF canopy 
area. Pioneer Aerospace (developers of the MPF 
parachute system) measured the area of existing MPF 
qualification parachutes. The flight and qualification 
parachutes were manufactured at the same time from the 
same lot of material. Hence, the final dimensions should 
be very similar. From these measurements, a mean 
value of the MPF parachute area and its variation were 
determined. This variation in parachute area was taken 
into account as part of the Monte Carlo methodology in 
determining the range on the MPF parachute Cd value. 

The values for the Backshell C D and Lander C D 
were also updated based on wind tunnel tests preformed 
at the ViGYAN low speed facility in Hampton, 
Virigina. 6 The tests improved the knowledge in the C D 
values of both configurations. 

Results 

Based on the reexamination of all the parameters 
described, the MPF parachute C D value is determined 
using a Monte Carlo analysis. The list of parameters 
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included in the analysis is given in Table 1, along with 
their mean value and uncertainty range. 

Using the specified distributions, 1000 random 
values were statistically sampled for each parameter 
within its uncertainty range. These values for all the 
parameters were then substituted in Eqn. (4) to 
determine an estimate for the MPF parachute C D . 

The resulting distribution in atmospheric density 
(using the procedure outlined in the Appendix) arising 
from the variation in atmospheric temperature and 
pressure is shown in Fig. 4. A variation of 
approximately +6% is observed about the mean density. 
The resulting distribution in the value of the 


MPF parachute C D arising from the variation in all the 
parameters is shown in Fig. 5. A mean value of 0.41 is 
calculated for the MPF parachute with a 3-0 uncertainty 
of +0.05 (±12%). A majority of this dispersion in the 
MPF parachute Cd value is due to the uncertainty in the 
atmospheric density. The actual dispersion in the Cd of 
a parachute is likely to be smaller than this range. 
However, due to the uncertainty in the atmospheric 
properties, a smaller range in the MPF parachute C D 
value is not possible. The minimum value obtained for 
the MPF parachute C D is 0.36. Note, in obtaining this 
revised prediction, conservative assumptions were 
utilized where appropriate. 


Table 1. Parameters included in Monte Carlo analysis. 


Parameter 

Mean 

Distribution 

Uncertainty Range 

m, kg 

520.9 

- 

- 

g, m/s 2 

3.7245 

- 

- 

Ab/si 

5.39 

- 

- 

^LaiD 

1.76 

- 

- 

Apar, m 2 

127.6 

Gaussian 

+5% (3-o) 

Vol, m 3 

135 

uniform 

+20% 

^ D B/S 

1.33 

uniform 

+5% 

C °Lan 

1.072 

uniform 

+5% 

Temperature, K 

221 

uniform 

+9 

Surface Pressure, mbar 

6.76 

uniform 

+0.15 

v, m/s 

65.5 

Gaussian 

+ 1.8 (3-o) 




Figure 4. Distribution in the Mean Atmospheric 
Density 


Figure 5. Distribution in the MPF Parachute C D value 
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Conclusions 

Based on the assumptions presented and 
reexamination of all the parameters (additional 
filtering of MPF radar altimeter data, revised 
atmospheric property estimates, and updated 
parachute area and Backshell and Lander Cd values), 
a revised flight reconstructed mean value for the 
MPF parachute Cd of 0.41 is obtained with a 3-0 
uncertainty of +0.05 (+12%). A majority of this 
variation in the MPF parachute C D value is due to the 
uncertainty in the atmospheric density. In obtaining 
this revised prediction, conservative assumptions 
were utilized where appropriate. 

The largest contributor to uncertainty in the MPF 
parachute Cd is uncertainty in density and the 
terminal velocity. Reexamination of the radar 
altimeter flight data allows for a better estimate on 
the terminal velocity and its uncertainty. In addition, 
re-examination of the atmospheric properties, 
supplemented by additional data, allows for a better 
estimate on the density and its uncertainty. 
Furthermore, the use of the Monte Carlo technique 
statistically allows for incorporation of uncertainties 
on all the parameters to better assess the variation in 
the MPF parachute C D value. 
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Appendix: Density Estimation Procedure 

In order to revise the Pathfinder parachute C D 
calculation during descent, an accurate prediction of 
the density of the atmosphere at some level near the 
surface is necessary, where the Pathfinder Lander is 
falling close to its terminal velocity. Density can be 
obtained from direct measurements, or can be 
inferred from measurements of pressure and 
temperature using Eqn. (A-l). 

The atmospheric density profile was derived 
from Pathfinder accelerometer measurements down 
to an altitude of just below 9 km, where the parachute 
deployed. Accurate measurements of surface pressure 
and temperature at 1 m were made immediately after 
landing, and at the same local time on subsequent 
days at the landing site. Unfortunately, no high 
quality temperature and pressure measurements were 
taken during parachute descent. 

In the original approach, an attempt was made to 
calculate the density at 300 m (altitude prior to airbag 
inflation) by extrapolating temperature and pressure 
from the known conditions at the surface. Pressure 
can be extrapolated fairly accurately using the 
hydrostatic equation (Eqs. (A-2) & (A-3)). To 
extrapolate temperature, results were used from a 1-D 
boundary layer model calculated for Pathfinder 
conditions, and which were shifted by a few degrees 
to agree with the 201.5 K Pathfinder measurement at 
1 m. 

Models show that for the Pathfinder landing 
conditions (Ls = 143°, LTST = 03:00, Lat = 19.5° N, 
Lon = 33.5° W), atmospheric temperature should 
increase from a minimum at the surface to a 
maximum at roughly 1 km, before falling with 
increasing altitude. There is considerable uncertainty 
in the width of the inversion layer below 1 km, so 
that the temperature extrapolation to 300 m might be 
subject to significant error. 
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MGS obtained radio occultation temperature 
profiles near the Pathfinder season, local time, and 
location, but for the following Mars year. These 
profiles have a vertical resolution of 0.5- 1.0 km, are 
accurate in the lower atmosphere, but do not extend 
completely to the surface. Four profiles selected 
within the bin Ls = 139.9°- 141.4°, LTST = 04:11, 
Lat = 16.5°- 20.9° N, Lon = 26.6°- 42.8° W, are 
consistent with each other to a few K below 10 km 
and show near surface temperature peaks of 220 K 
near 1 km altitude. 

The MGS radio-occultation data suggested a new 
approach. Densities calculated near 1 km (the 
atmospheric temperature peak) are less sensitive to 
vertical gradients, and can be based on reliable radio- 
occultation temperature measurements. The 
recommended density calculation utilized in this 
analysis is summarized as follows: 

1. Derive p(z) from P(z) and T(z) at 1 km, using 
Eqn. (A-l) . 

2. Based on the occultation data, T(z) = 221 K, 
+3K(l-o) at lkm. 

3. P(z) at 1 km is derived from Eqs. (A-2) & 
(A-3). 

Given: Ps = 6.76 +0.05 mbar (1-a) from Pathfinder 
(Eqn. (A-2)). 

T = 215 Kin Eqn. (A-3). 

Note: T is the mean temperature from 0-1 km, which 
is probably weighted more towards the 1 km than the 
1 m temperature because the surface temperature 
inversion is significantly narrower than 1 km. The 
consequences of this are not very significant in the 
pressure calculation. 

From the ideal gas law, the density at altitude z 
km in Mars’ atmosphere is given by the expression: 

p(z) = [P(z)/P0] * [T0/T(z)] * [M/Vm] (A-l) 
where 

p(z) is the density at altitude z in gm/cm 3 
P(z) is the pressure at altitude z in mbars 
T(z) is the temperature at altitude z in Kelvin 
P0 = 1013.24 mbar (Standard Pressure) 

TO = 273.15 K (Standard Temperature) 

M = 43.5 gm (Mean atmospheric molar weight) 
Vm = 2.241e+4 cm 3 (Molar volume at STP) 

Given P(z) and T(z) the density can be calculated 
Eqn. (A-l). T(z) at the level of interest must be 


specified. Provided z < 2 km, P(z) can be 
extrapolated from the surface using the expression: 

P(z) = Ps * Exp(-z/H) (A-2) 

where 

Ps is Surface Pressure 

H is the mean atmospheric pressure scale height 
in km 

H is approximated by the expression: 

H = RT/(Mg) (A-3) 

where 

R = 8.314 J/K/Mole (The gas constant) 

T = Mean temperature from 0-z km (need not be 
very accurate) 

M = 43.5 gm (Mean molar weight as in Eqn. (A-l)) 
g = 3.717 m/s 2 (Acceleration due to gravity at 
surface) 
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